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Abstract 

By imposing axial and vector Ward identities for flavour-non-singlet 
t^. . currents, we estimate in the quenched approximation the non-per- 

turbative values of combinations of improvement coefficients, which 
appear in the expansion around the massless case of the renormaliza- 
t*"*- ■ tion constants of axial, pseudoscalar, vector, scalar non-singlet cur- 

rents and of the renormalized mass. These coefficients are relevant 
for the completion of the improvement programme to 0(a) of such 
operators. The simulations are performed with a clover Wilson action 
non-perturbatively improved. 
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1 Introduction 



The programme of the improvement jl] of the Wilson action has been actively 
developed at the non-perturbative level over the last years 0, |, ||. At first, 
in the framework of the Schrodinger functional it was possible to determine 
non perturbatively the dependence upon the bare coupling constant of the 
coefficient csw @ of the clover term in the improved action || . 

Besides improving the action, the programme includes the improvement 
of the operators appearing in the correlation functions related to phenomeno- 
logical interesting quantities such as pseudoscalar meson decay constants and 
the matrix element of the four-fermion operators of the weak effective Hamil- 
tonian. 

In general the operator improvement consists of two parts: the mix- 
ing with higher-dimensional operators with the same quantum numbers (in 
the literature the mixing coefficients are called c) and the multiplication by 
a suitable renormalization constant. The ultraviolet-finite renormalization 
constants can be expanded around the massless case: 

Z (m ^ 0) = Z ■ (I + b ma + ...), (1) 

where we have omitted corrections due to lattice artefacts of order a 2 and 
higher. 

Some of these quantities have been calculated at the perturbative level 
for axial (A), vector (V), pseudoscalar (P) and scalar (S) currents as well 
as for the renormalized mass (m) 0, |8|, ||: non perturbative estimates are 
available for Za, Z v , by |J and for ca @ and cy . 

In this letter, we present a non-perturbative determination of the quan- 
tities: 

b A - bp, by - b s , b m (2) 

and of the ratios: 

Z m Zp/ZA, Z m Z s jZy (3) 

from a set of axial and vector Ward identities. 



2 The method 

The extraction of the b coefficients and of the Z ratios is based on the fol- 
lowing Ward identities: 

d^{x)n\0)) = 2m jk (P(x)^( ))+O(a 2 ) 
9 M (Vj(x)fi+(0)) = Am jk (S(x)fi t (0))+O(a 2 ), (4) 

which, after x integration, become: 

dt(Ai{t)Q\0)) = 2m jk (P(t)Q\0))+O(a 2 ) 
dtiVjWrtiQ)) = Am jk (S(t)Q\0))+O(a 2 ). (5) 
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The suffix / for the axial current indicates that the current is improved by 
the appropriate mixing with the pseudoscalar density multiplied by the coef- 
ficient ca- The value for ca is taken from its non-perturbative determination 
in ref.pl. For the vector current, the contribution of the mixing with the 



tensor current (see ref. |10| for non-perturbative determination of the mixing 
coefficient cy) vanishes because of the antisymmetry of tensor indices. 

The indices j, k refer to the flavour content of the bilinear operators, 
which can be written as 

o(t) = £o(x) =Y,h(x)roMx)- ( 6 ) 

X X 

Equations hold for any operator Q at any time different from t, 
reflecting the fact that the W.I. are identities among operators. For the axial 
W.I. we use Q = P(0) and Q = Aq(0) while for the vector W.I. we use 
Q = S(0) and Q = Vo(0); in both cases the latter operator leads to much 
noisier results. 

The quantities rrij k and Arrij k are related respectively to the average 
and the difference of the renormalized mass, according to 

R Zp 1 + bpfnZa 
m jk = m jk —-—— = 
J Z A 1 + b A m q a 



Arn jfc - Am jk — (( 
J Z v 1 + bym q a 



m fk = \i m f + m%), Amf k = mf - mf , 



where m q a is the average of the bare masses j and k: 



1, 

m q a = -{m q] a + m qk a) 



with k c the critical value of the Wilson hopping parameter. The parameter 
k c is determined from the chiral extrapolation of the mass defined through 
the axial Ward identities themselves. 

The renormalization constants can be determined by replacing the "un- 
renormalized current masses" mj k and Arrij k with the renormalized ones 
through the above equation, and then the renormalized masses in terms of 
the bare masses: 

m R = Z m m q (l + b m m q a). (9) 
Indeed, by including the lattice artefacts up to 0(a): 

m jk a = Zp J* m ffn^a - (b A - b P )(jn^af + b m (m q a) 2 ) 

Zj A V ' 

Z Z 

Am jk a = m (m qj a - m qk a) (1 + (2b m - (b v - 6 5 ))m^a) (10) 
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From a fit of the bare mass dependence of these results we can determine 
non-perturbatively the combinations in eqs. (0,0). 

Within the approximation used in the above formulae, we can derive 
an expression for (pA — bp): 

(rrijja + m kk a) Z P Z m 1 2 2 

m jfc a — -{b A - bp){m qi - m qj ) a (11J 

which can be used to determine a value for the combination independently 
from the knowledge of the critical value of k. We have compared such a 
determination of b A — bp with the one coming from a global fit to expression 
eq. QlOD and used it as a sensitive check of our estimate of k c . 

The presence of flavour-non-diagonal currents is important for the fit 
of the axial W.I. and essential for the vector ones. We want to point out that 
rrijka depends either upon (m^a) 2 or upon (m q a) 2 , which are different vari- 
ables if the flavours are not degenerate, making then possible to disentangle 
the coefficient b m from the combination 6a — bp. 

Our fits are to the dependence upon the quark mass, and order a 2 
corrections linear in the quark mass can in general fake the extraction of 
the coefficients. The simple lattice discretization of the time derivative 
j-(f(t+a)-f(t-a)) has an error of the order of f^a 2 , which becomes f^M 2 a 2 
when a single state of mass M dominates the correlation function /. In the 
pseudoscalar channel, chiral symmetry makes this term proportional to the 
quark mass. In other channels, the mass M acquires a non-zero value when 
the quark masses vanish. The quantity (Ma) 2 still contains a term that is 
linear in the quark mass but not in the lattice spacing. In order to minimize 
such effects we have used the lattice discretization of the time derivative ad t 
correct up to term f^a 5 . While for the pseudoscalar case with spontaneous 
symmetry breaking the improvement to the fifth order of the derivative re- 
moves these fake linear terms in the quark mass, for the vector current a 
linear term survives at any finite order n, although suppressed by a coeffi- 
cient of order l/(n — 1)!. We have checked that improving the derivative to 
the next order does not change our results beyond their accuracy. 

We cannot exclude the presence of other lattice artefacts in the ratio 
of matrix elements formally of order a 2 but linear in the quark mass. We 
have checked in the case of the axial current that the results are stable with 
respect to the choice of the operators and we interpret this as a sign for the 
absence of large extra terms of order mAa 2 . 

This method of computing the combinations in eqs. (0H) is valid in 
the quenched approximation. The presence of dynamical quark loops would 
introduce an additional sea quark mass dependence, which would involve 
flavours different from those in the currents. 



3 The numerical results 

We have performed several simulations at different values of f3, in order to 
derive the coupling constant dependence of the quantities in eqs. (0,0). The 
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Table 1: The values of k used in our simulations at various (3 



values of (3 used in the simulations and the corresponding volumes are col- 
lected in Table |I], with a list of the values of k and of our best estimate of 
the critical k obtained from the W.I. themselves. The values of well 
compatible within errors with those of ref. ||. The variation of our results 
under a change of k c within the quoted error is smaller than the accuracy by 
which we can extract the non-perturbative quantities from our fit. 

Simulations are performed with an updating sequence made by a stan- 
dard heat-bath followed by 3 over-relaxation steps. The improved fermion 
propagator is calculated every 1000 gauge update using a stabilized bicon- 
jugate algorithm. For our runs we have used the 25 Gflops machine of the 
APE series, made of 512 nodes working in SIMD (Single Instruction Multiple 
Data) mode. 

Each propagator was summed over the space volume distributed to the 
single node (3x3x2) and stored on disk. The use of these propagators 
leads to correlation functions that contain the correct local-gauge-invariant 
terms and other non-local, gauge-non-invariant terms that go to zero after 
summing over the gauge configurations. We have explicitly checked that 
with our statistics the residual noise due to imperfect cancellation of the 
gauge-non-invariant terms is much below the statistical fluctuations. Storing 
fermion propagators allows for an off-line calculation of all flavour-non-singlet 
correlations. 

The W.I. are satisfied separately at each time: after some initial time, 
up to which higher-order lattice artefacts still dominate, rrijk and Arrijk show 
a plateau. At (3 — 6.2 we run two temporal extensions (32 and 48) in order to 
monitor the stability of the plateau. We have used two methods of analysis: 
either we first average the result over the time interval of the plateau and 



4 



r3 t 

Li 1 


10 4o 


10 t)Z 


10 OZ 


10 OZ 


a 
P 


P. 9 
0.Z 


O.o 


o.U 


1 9 n 
1Z.U 


7^ confs 


ou 


ou 


oU 


oU 


ZtuZp/Za 


1.09(1) 


1.08(1) 


1.08(1) 


1.060(6) 




1 9/Ll'9"\ 

1 Z 1 


l.Zll^lj 


1 . 14Z 14 1 


i.uou^o j 


b A - bp 


0.15(2) 


0.10(2) 


0.06(2) 


0.04(2) 


b 


-0.62(3) 


-0.58(3) 


-0.57(3) 


-0.53(2) 


b m - (by - b S )/2 


-0.69(4) 


-0.63(4) 


-0.54(3) 


-0.52(2) 



Table 2: The results of our calculations 



then perform a fit, or we first perform a fit at each time value inside the 
plateau and then average the results of the fit at different times. The two 
procedures give very similar results. 

The choice of the spectator operator Q affects the statistical error of the 
final results. We have found that the pseudoscalar and the scalar densities 
give the best results for the axial and the vector case respectively. 

We perform a fit of 2rrijk with the function (here all masses are in lattice 
spacing units): 

ai(m qi +m qj ) +a 2 (m 2 q . +m 2 q , ) +a 3 (m qi +m qj ) 2 +a 4 {m\. +m q , ) +a 5 m qt m qj (m qi +m q . ) 

The first three coefficients of the fit are related to the renormalization con- 
stants as follows: ai = Z m Zp/ZA', a 2 /a\ = b m ; a 3 /ai = —(bA — bp)/2. 
The last two coefficients in the fit can be introduced to parametrize order- 
a 2 corrections compatible with the flavour exchange symmetry of the Ward 
identity. 

For the quantity Am^ we perform the fit with: 

v x {m qi - m qj ) + v 2 (m 2 g . - m 2 q .) + v 3 (m q . - m 3 q .) + v 4 jn qi m q .{m qi - m q .) 

where V\ = Z m Zs/Z v ; v 2 fv\ = b m — (by — bp)/2. As before, the extra 
coefficients t>3 and V4 parametrize the possible order-a 2 corrections. 

Our results normally refer to the fit with three parameters for the axial 
case and two for the vector case. Increasing the number of parameters in 
general does not improve the value of x 2 much while it considerably increases 
the error. The results are compatible with the lower parameter fit, with the 
exception of the determination of by, which comes systematically higher with 
the four-parameter fit. We have included this effect in the corresponding 
error. 

Table ^| contains the main results, the values for the various renormal- 
ization parameters at different (3 and volumes. The f3 = 6.2 results on the 
smaller temporal extension are fully compatible. 

With the values of the fit we can check if the renormalized W.I. depend 
only upon the sum (for the axial) or the difference (for the vector) of the 
renormalized masses. 
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Figure 1: The non-perturbative result for b A — b P . The perturbative result 
of 0(g 2 ) is negligible on this scale. 
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Figure 2: The non-perturbative estimate of b m is compared with the pertur- 
bative result. 
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Figure 3: The non-perturbative result for by is compared, after using 
Liischer's relation, with the one of ref. [6] (dotted curve) and with the per- 
turbative result (dashed curve). 
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Figure 4: The non-perturbative result for Z v is compared, after using 
Liischer's relation, with the one of ref. [6]. 
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The renormalized masses and currents manage to bring on the same 
straight line points that appear misaligned and on a curved line for the bare 
quantities. For large values of the masses, higher-order terms enter the game 
and produce again a misalignment of the corresponding points. 

The fits in general do not include all k values; we exclude the heavier 
masses until the stability of the results is reached. 

Our results, when compared with those available from perturbation 
theory, show that higher-order terms seem to dominate for the differences 
bA — bp, which is very small at one-loop order (see fig. [1]), while for b m 
the presence of sizeable terms of order g 2 makes the effect of g 4 terms less 
prominent. Indeed, our results are not far from lowest-order perturbation 
theory for this case (see fig. 0). 

For by — bs + 2b m , there is an argument due to Martin Luscher [11] relating 
b m with bs and Zs with Z m in the quenched approximation: 2b m — bs = 
and ZsZ m = 1, which implies that from the W.I. for the vector current 
we actually obtain by and Zy. The comparison with those of ref. |J is 
shown in figs. || and |]: while for Zy there is a perfect agreement, for by we 
are generally closer to the perturbative result. Our large errors are mainly 
systematic and reflect the instability of a four-parameter fit. Running at 
lower quark masses could reduce the discrepancy which might also be due to 
residual order-a 2 lattice artefacts. 



The use of axial and vector Ward identities with flavour-non-singlet 
currents allows the determination in the quenched approximation of various 
non-perturbative renormalization constants. The calculation that we have 
presented could be refined by using the Schrodinger functional method which 
would allow a safe investigation of the very low quark mass region. 
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